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Abstract 

A scheme to calculate the electronic structure of systems having a spiral magnetic structure is 
presented. The approach is based on the KKR (Korringa-Kohn-Rostoker) Green's function formal- 
ism which allows in combination with CPA (Coherent Potential Approximation) alloy theory to 
deal with chemically disordered materials. It is applied to the magnetic random alloys Fe x Nii_ x , 
Fe^Coi-z and ¥e x yha.i-. x . For these systems the stability of their magnetic structure was ana- 
lyzed. For Fe^Nii-a; the spin stiffness for was determined as a function of concentration that was 
found in satisfying agreement with experiment. Performing spin spiral calculations the longitudinal 
momentum-dependent magnetic susceptibility was calculated for pure elemental systems (Cr, Ni) 
being in non-magnetic state as well as for random alloys (Ag x Pti_ x ). The obtained susceptibility 
was used to analyze the stability of the paramagnetic state of these systems. 

PACS numbers: 71.15.-m,71.55.Ak, 75.30.Ds 
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I. INTRODUCTION 



The use of symmetry properties of solids for calculations of their electronic structure is a 
very efficient way to reduce the computational effort required for the solution of the problem. 
In particular, the single-particle electronic states of paramagnetic or collinear magnetic infi- 
nite solids can be effectively found by solving the corresponding Kohn-Sham-Dirac equation 
making use of the Bloch theorem. Dealing with systems exhibiting non-collinear magnetic 
structure, the electronic structure problem becomes much more complicated because of bro- 
ken symmetry (in general, both translational and rotational), leading to an increase of the 
unit cell of a system and a corresponding increase of the required computational effort. 

Sandratskii introduced an approach that allows to calculate the electronic structure of 
systems with spiral magnetic structures in an efficient way . This approach is based 
on the symmetry properties of spin spiral structures as investigated by Brinkman and Elliot 
3, 4 1 and Herring j^j and allows to deal with long-period non-collinear magnetic structures 
avoiding the use of big unit cells in electronic structure calculations 6]. This makes it an 
efficient tool for the analysis of the stability of various non-collinear magnetic structures 
with different translation period, as for example demonstrated by Mryasov et al. |7| for the 
investigation of the magnetic structure of fee Fe. 

In the case of systems with a collinear magnetic structure as a ground state spin spirals 
can be treated as transverse spin fluctuations in the adiabatic approximation. The energy 
dispersion of such fluctuations AE(q) gives access to the spin stiffness and exchange coupling 
constants of a system and in this way to the spin excitation spectrum as well as finite 
temperature magnetism 8MlO| . An important feature of spin spiral calculations is that they 
account for longitudinal fluctuations of the magnetic moment. This leads to more reliable 
results for AE(q) compared to those obtained using the non-self consistent force-theorem 
approach. 

As was pointed out by Sandratskii and Kiibler [11] the technique for spin spiral calcu- 
lations can be used for calculations of the static (u> = 0) momentum-resolved longitudinal 
magnetic susceptibility. Untill now only few corresponding ab-initio calculations have been 
presented in the literature. In most cases the static g-dependent magnetic susceptibility was 
calculated using perturbation theory [ijl or performing super-cell calculations jl3]. The 
spin spiral method, on the other hand, allows to perform self-consistent calculations of the 
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magnetic susceptibility avoiding the super-cell concept 11]. 

All spin spiral calculations were done so far using the ASW [H 
band structure methods. These methods use Bloch-function basis sets to represent the 
solution of the Kohn-Sham equation and for that reason are restricted to ordered materials 
concerning then application. Use of multiple scattering theory in combination with CPA 
(Coherent Potential Approximation) alloy theory, on the other hand, substantially extends 
the variety of materials which can be investigated by giving access to systems without 
chemical order. Here we present the implementation of spin spiral approach within the 
Korringa-Kohn-Rostoker (KKR) Green's function band structure method We will 

show results of calculations for different systems focusing on disordered alloys. 



II. THEORETICAL BACKGROUND 

When dealing with the electronic structure of solid state systems having a spiral mag- 
netic structure rotations can be applied independently to the spin and spatial parts of the 
electronic wave function if spin-orbit coupling (SOC) is neglected. Using a spin-diagonal 
form of the exchange-correlation potential in the local frame of reference of an atom site, 
the Kohn-Sham equation for the spinor wavefuction ip(r) can be written in the form: 
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ip(r) = Eip(r) . (1) 



Here r v denotes a position of an atom in a unit cell, R n is a Bravais lattice vector and U nv 
is a spin transformation matrix that connects the global frame of reference of the crystal to 
the local frame of the atom site at r v + R n that has its magnetic moment tilted away from 
the global z-direction. The transformation U nu is characterized by the Euler angles 9 njU and 
<f) nM as is is shown in Fig. [1] for the case of a spin spiral. 

As was shown by Sandratskii, considering a spin spiral structure, Eq. ([1]) can be easily 
dealt with using the properties of spin space groups (SSG) allowing independent trans- 
formations within the spin and space sub-spaces. The spin spirals characterized by the wave 
vector q angles 9 U and (p u are represented by the expression: 

m v n = m v [cos(<f • R n + (p u )sm9 u , sm(q ■ R n + (p u )sm9 u , cos^] 
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FIG. 1: Geometry of a spin spiral with the wave vector q along the z-direction. 



defining the spin direction at every site (n, v) of the lattice with m u the magnitude of the 
magnetic moment on site v within the unit cell. Assuming a collinear alignment of the spin 
density within the atomic cell at (n, v), it's natural to use a local frame of reference with its 
z-axis oriented along fh v n . The corresponding transformation matrices U w occurring in Eq. 
([1]) can be written as a product of two independent rotation matrices U^ v = U n (9 UA 4>u, q) = 
U V (9 U , 4>v)U^ n where the matrix depends only on the translation vector -R n :[l|, |2| 
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Instead of solving the Kohn-Sham Eq. ([T]) for the eigen functions and values the electronic 
structure can be represented in terms of the corresponding Green's function. Within multiple 
scattering theory the Green's function is represented in real space by the scattering path 
operator f nn together with the regular Z% (r,E) and irregular J^if^E) solutions of the 
single-site Kohn-Sham equation referring to the local frame of reference: 



G + (r,r',E) = £ Z^f, E)r^' A2 (E)Zl*(r E) 

AiA 2 

-Y,[zt(r,E)J^(f',E)Q(r'-r) 

Ai 

Jl(r,E)Zl*(r>,E)e(r-r') 



(3) 
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The scattering path operator is defined by its equation of motion: 

•nvn'v' 4/nv x i v ' s~inqmu, rniin' V / a\ 

Z = t nun > u > + t 7 ^ G Z ; ( 4 ) 

kp 

where t nu and G nqmfl are the single-site t-matrix and free-electron propagator, respectively, 
that are all expressed with respect to a common global frame of reference. Eq. (J3j) has the 
formal solution 

Z=\f 1 -^ 1 . (5) 

In Eqs. (H]) and (jHJ) the underline indicates matrices in the (I, m) representation while dou- 
ble underline indicates super-matrices including the site index. In the case of a collinear 
magnetic structure the local and global frames of reference coincide. This implies that Eq. 

gives immediately the solution with respect to the local frame of reference. For infinite 
systems having a regular periodic lattice a solution to Eq. (jlj) can be obtained by Fourier 
transformation instead of using the real space expression given in Eq. (J5j) . 

For non-collinear magnetic solids with a periodic lattice structure one can solve Eq. (jlj) 
as for collinear systems but using an extended super-cell. The size of the corresponding 
super-cell is determined by the period of magnetic structure. All atoms within the cell are 
in general inequivalent and have their own local frame of reference. Therefore super-cell 
calculations can be rather time consuming in particular for magnetic structures having a 
long period. 

However, as pointed by various authors js-5] use of symmetry allows to simplify the 
problem substantially. Spiral magnetic structures transform according to the group of gen- 
eralized translations that are characterized by the wave vector q and represented by the 
matrices (Eq. ([2])). This implies in particular that the matrices allow to express 
the single-site t-matrix t nv at site {nu) to that at site (Ou). This symmetry property allows 
to write the scattering path operator referring the global frame of reference as follows: 

•nvn'v 1 4.nv r , iTi!>\ v ' /~<nv mil rutin 1 v 1 

T = t Onvn<v>+t 2^ 

mfj, 

= Uffi w Ui5 ni/n , u , + U^t w UiY^G nvm ^r m ^' . (6) 
This allows to find the scattering path operator and from this the Green's function in the 
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local frame of reference of the each atom, solving the equation 
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where the tilde indicates matrices which refer to the local frame of reference. 

In the last line of Eq. (JSJ) use has been made that the single-site t-matrices t do not 
depend on the lattice index n but only on the site index v in the unit cell. As a consequence, 
the multiple scattering problem can be solved as for the case of collinear magnetic structures 
by Fourier transformation of the equation of motion for the scattering path operator. This 
leads to its representation in reciprocal space according to: 
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The structural Green's function referring to the local frame of reference can be determined 
as follows: 
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(10) 



Here G vv (k) is a structural Green's function for one spin channel represented in the global 
frame of reference. 

The charge distribution within the central unit cell n = is determined by the cell- 
diagonal scattering path operator r 00 which can be found by the Brillouin zone integral 
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where Uq is the transformation matrix diagonalising the potentials as well as i-matrices with 
respect to spin within the central unit cell. 

To perform calculations for disordered alloys the CPA (Coherent Potential Approxima- 
tion) alloy theory [l?} [isl is used. In the case of a spin spiral system the CPA medium is 
represented in the global frame of reference by the effective single-site scattering matrix t 
and the scattering path operator obtained from the expression: 
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r 00 ' c (£) = — - / d k (t (E)) — G(k, E) . (12) 
^bz Jn BZ L J 

The corresponding element projected scattering path operators are obtained from these via: 
with 

R a = [1+ [(f)- 1 - if)' 1 ] z 00 ' c ] 1 • (14) 

The approach developed for calculations of non-collinear spin spiral structures can be used 
for investigations on the longitudinal magnetic susceptibility as a function of the wave vector 



q ll( . This approach allows to avoid the use of perturbation theory and can be applied to 
magnetic as well as non-magnetic systems. In the following we focus on materials in their 
non-magnetic state which may exhibit paramagnetism (Ag x Pti_ x ), ferromagnetism (Ni) or 
antiferromagnetism (Cr) in their ground state. For this purpose we specify a spiral external 
magnetic field to be perpendicular to the direction of the wave vector q (i.e. 9 = 90°): 

h(f) = h [cos(q- Rn), sin(<f • R n ), 0] . 

In this case is the potential energy term in the Kohn-Sham equation (see Eq. (flj)) is given 
by: 

\ V nv {r) + B%* + A^ c (f) 1 



nv \ u v i y'v j 



(15) 

A self-consistent calculation based on Eq. (|15p gives the spin magnetic moment induced by 
the external magnetic field. The <f-dependent external magnetic field should be taken small 
enough to be considered as a perturbation. In this case, assuming a linear response to be 
the leading term of the response function the corresponding magnetic susceptibility can be 
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derived from the expression 



Xkl = ^ ■ (16) 
no 



Suppressing the spin-dependent part of the exchange-correlation potential (AB^ c (r) = 0), 
one can calculate the unenhanced spin susceptibility x°(o)- Otherwise, Eq. (JTBjl gives the 
enhanced longitudinal magnetic susceptibility x(o)i represented in linear response theory for 
uniform system by the expression 

x(q) = = pPrr= , (17) 

with I(q) the exchange integral responsible for the enhancement of the magnetic suscep- 



tibility (see, for example, 10|, |l9|). In case of a paramagnetic ground state the magnetic 



susceptibility x{q) is positive for all values of q. For other cases the denominator in the Eq. 
( IT7|l may become zero or even negative. This singular behavior of the susceptibility obvi- 
ously indicates an instability of the paramagnetic state towards a transition to spontaneous 
formation of ferro- or anti-ferronagnetic order. 



III. RESULTS 



A. Spin spiral structure in alloys 

In the following several applications of the scheme introduced above are presented that 
focus on disordered alloys to demonstrate the flexibility of the multiple scattering formalism 
when dealing with spin spiral systems. Corresponding calculations have been performed for 
alloys having fee (Fe^Nii-^ and Fe^Mni-^) and bec (Fe^Co!^) lattice structure. 

For all calculations the angle 9 has been chosen to be 90°. For this spin geometry, the 
spin spiral with q = ^(0,0, 1) corresponds to a spin configuration where the first neighbor 
atoms in the (001) direction have an anti-parallel (AFM) spin alignment, while q = implies 
a parallel (FM) orientation. 



1. The disordered alloy system Fe x Ni\_ x 

Fig. [2k. shows the energy of the disordered Fe^M^a. alloy system with a spin spiral 
structure as a function of the wave vector q. For all concentrations the minimum of the energy 
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FIG. 2: a) The energy of spin spiral magnetic structure in Fe x Nii_ x alloys and b) local magnetic 
moments on Ni atoms as a function of the wave vector q= ~(0, 0, q z ). 



is found for q = 0, implying that the ferromagnetic structure is more stable configuration 
than non-collinear structures characterized by wave vectors along the (001) direction. 

As can be seen from Fig. Eb, the local magnetic moment of Ni in pure Ni drops down to 
m = at the wave vector q c « ^(0, 0, 0.5) and the system becomes paramagnetic. In terms 
of the Stoner theory of ferromagnetism (see, e.g. Ref. ilQl ) this means that the criterium for 
the instability of the paramagnetic state is satisfied only for small wave vectors, while above 
q c the paramagnetic (PM) state should be the most stable state of the system. The criterion 
for the instability of the PM state will be discussed below in more detail. 

Adding only small amounts of Fe to Ni leads obviously to a nonzero magnetic moment per 
unit cell at all values of wave vector q. This is caused by the large magnetic moment of Fe 
which depends only slightly on the wave vector. Fig. [2b shows that the Ni magnetic moment 
in contrast to that of Fe, varies rather rapidly with increasing wave vector and changes sign 
at q ~ ^(0, 0, 0.6). This means that in the vicinity of the ground state of the alloys (q = 0) 
the magnetic moments of Fe and Ni atoms prefer to have parallel alignment, while close to 
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FIG. 3: Spin stiffness constant of Fe^Mi-^ alloys as a function of the concentration in comparison 
with experiment: Nakai [2(| (full circles), Hatherly et al. [21 1 (open circles), Rusov |23 | (full 
squares). 

q — ^(0, 0, 1) (AFM structure along (0,0,1) direction) the more favorable orientation of the 
Fe and Ni moments is anti-parallel (AP). Nevertheless, even for small Fe concentrations, the 
total magnetic moment is determined by the dominating moment of Fe. As a result, the 
alloy system exhibits effectively a ferromagnetic behavior for all wave vectors, as one can 
see in Fig. [2j 

The energy difference between the spin spiral states with q = and q = ^(0, 0, 1) remains 
almost unchanged up to the Fe concentration XFe ~ 0.6, and changes nearly by 20 % when 
approaching xp e ~ 0.65. On the other hand, the spin-stiffness constant deduced from 
the energy dispersion curves decreases continuously with the increase of Fe content as can 
be seen from Fig. [3j This figure also shows that the calculations reproduce the available 
experimental data for the spin-stiffness constant fairly well, although they seem to be slightly 
to high. This difference can be partially attributed to the conditions of the experiment as 
e.g. polycrystallinity of the samples and a finite temperature. 



2. The disordered alloy Feo^Coo^ 

The change of sign of the magnetic moment observed for Fe^Mi-^ alloys for one of the 
alloy components upon variation of the wave vector becomes even more pronounced in bcc 
Fe .5Co .5 and fee Fe .5Mn . 5 alloys. Disordered bcc Feo.sCoo.s has a ferromagnetic ground 
state. The spin spiral energy shown in Fig. increases with wave vector confirming the 
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FIG. 4: a) the energy of spin spiral magnetic structure in Feo.sCo-o.s calculated for the wave vector 
q = ^(0, 0, q z ) along [001] direction; b) local magnetic moments on Fe and Co atoms separately, as 
a function of wave vector of spin spirals. 

stability of the FM state. As can be seen from Fig. around q = the individual Fe 
and Co moments are aligned parallel with respect to each other. However, after crossing 
q c ps ^(0,0,0.7), the total magnetic moment jumps from m = 1.46/ie to m = 0.47/x.b due 
to a change of the sign of the Co magnetic moment with respect to that of the dominating 
Fe moment. As Fig. H)d shows, the dispersion of the spin spiral energy for the anti-parallel 
configuration gets very weak up to q ~ ^(0,0,1). To estimate the energy of the spin 
spirals for the non-equilibrium configurations, i.e. anti-parallel for q c < ^(0,0,0.7) and 
parallel for q c > ^(0, 0, 0.7), respectively, frozen potential calculations have been performed. 
The corresponding results are represented in Fig. 0b by dashed and dashed-dotted lines. 
Obviously, these results augment the two stable branches fairly will. 



3. The disordered alloy Feo.^MriQ^ 

Fig. shows the results of spin spiral calculations for Feo.5Mn .5 having a non-collinear 



magnetic structure as a ground state [23H25I] . As is seen in the energy dispersion curve, Fig. 
[5^, the system exhibits the behavior of a FM system for wave vectors q up to q c ps ^(0, 0, 0.6). 
In this wave vector region the alloy has a small average magnetic moment formed by two 
anti-parallel aligned magnetic moments of Fe and Mn occupying randomly the site (see Fig. 
Eb). 

At q c ~ ^(0,0,0.6) the energy of a spin spiral reaches its maximum and the following 
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FIG. 5: a) The energy of spin spiral magnetic structure in Fe^M^-^ alloys calculated for the 
wave vector q = ^-(0,0, q z ) along [001] direction; b) local magnetic moments on Fe and Mn atoms 
separately, as a function of wave vector of spin spirals. 

increase of the wave vector is accompanied by a decrease in energy and an increase of the 
average magnetic moment. At q = (0, 0, -) the spin spiral magnetic structure reaches its 
energy minimum, which is about 50 meV lower than the energy of the FM state, with a 
parallel alignment of the magnetic moments of the alloy components. 

Similar to Feo.sCoo.s, these two minima of the energy - around q = and around q = 
(0,0,^) - are formed by two crossing branches of the spin spiral dispersion relation: one 
corresponds to an anti-parallel alignment of the Fe and Mn magnetic moments (around the 
q — 0) and another to their parallel alignment (around q = (0, 0, ~) ), which have a crossing 
point at q « * (0,0,0.6). 

Thus, from the analysis of the energetics of the spin spiral structures in Feo.5Mn .5, one 
can conclude that the system has in its magnetic ground state an anti-parallel alignment of 
the magnetic moments of first neighbors, no matter whether the neighboring atoms are Fe 
or Mn. 



B. Spin susceptibility 

In the present section we will discuss another application of the technique presented above. 
As was shown by Sandratskii and Kiibler spin spiral calculations can also be used to 
determine the longitudinal magnetic susceptibility x, both for magnetic and non-magnetic 
systems, as a function of the wave vector q. This approach allows in particular to avoid 
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FIG. 6: a) Wave-vector (q = f (0, 0, g z )) dependent magnetic susceptibility of the paramagnetic 
disordered Ag^Pti-^ for various concentrations, b) Comparison of the present results for the 
susceptibility for q = with the results of Deng et al. [fjj] obtained via linear response theory. 




the use of perturbation theory. Adding a Zeeman term to the Hamiltonian corresponding 
to a small external helical magnetic field allows to obtain the magnetic susceptibility from 
the the induced magnetic moments. For the present calculations a Zeeman splitting ho — 1 
meV has been used. 

The present work deals with non-magnetic systems, which have either a paramagnetic 
(AgPt), a ferromagnetic (Ni) or an anti/ferromagnetic (Cr) ground state. Dealing with 
magnetic systems being in an imposed paramagnetic state, their magnetic susceptibility 
gives information on an instability with respect to magnetic ordering. 

1. The paramagnetic disordered alloy Ag x Pt\^ x 

Fig. [6] shows the magnetic susceptibility of paramagnetic Ag^Pti-z alloys as a func- 
tion of the wave vector q for various concentrations. The spin susceptibility of the al- 
loys presented in Fig. [6] is composed by contributions from both components according to 
x(o) = x XA g (q,x) + (1 — %)XPt(<l, x). For all concentrations, the increase of the wave vector 
for helical magnetic field is accompanied by a decrease of the response functions, as it is 
usually found for paramagnetic systems. The main contribution to the spin susceptibility 
stems from the Pt atoms. As can be seen, increasing the Ag content leads to a decrease of 
the magnetic susceptibility for all values of wave vector. 

The present results for q = are compared with the total magnetic susceptibility obtained 
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via fully relativistic linear response calculations [19| . As one can see, the agreement of results 
obtained by the two rather different theoretical approaches is rather good. One reason for 
the observed deviations is the use of a finite value for the external magnetic field in the 
present calculations giving the magnetic susceptibility from the induced magnetic moment 
within the self-consistent calculations. Another reason is the neglect of spin-orbit coupling 
within the present calculations that usually reduces the spin susceptibility. Nevertheless, 
both approaches lead obviously to coherent results that are in rather satisfying agreement 



with experimental results 26] (full squares in Fig. |6b). Note however, that experimental 



results represent the total magnetic susceptibility including also the orbital contribution. 

2. Pure ferromagnetic fee Ni 

The calculations performed for ferromagnetic Ni in a paramagnetic state show a behavior 
for the magnetic susceptibility as a function of the wave vector that is rather different from 
that of systems with a paramagnetic ground state as for example Ag x Pti_ x alloys) (see 
Fig. [7]). The paramagnetic state of Ni was simulated using the disordered local moment 
(DLM) 23] method assuming equal concentration for atoms with opposite orientation of 
their magnetic moments. The magnetically disordered state of Ni is characterized by a 
vanishing local magnetic moment and therefore the DLM method allows us to force the 
local magnetic moment to be zero. Fig. [7Ji shows the results obtained for Ni with the 
experimental lattice parameter a = 6.65 a.u. At small values of the wave vector q the 
magnetic susceptibility is negative indicating an instability of the paramagnetic state. This 
is a result of the high density of states (DOS) of the 3d-electrons leading to a large value 
of the unenhanced magnetic susceptibility x°- Accordingly, for small (f-vectors the Stoner 
condition for a mag netic instability I(q)x°(o) > 1 ( E Q- [IB (see, e.g., |28|,|29|) is fulfilled. 

As one can see in Fig. [71 at the wave vector q « ^(0, 0, 0.55) (for which the denominator in 
Eq. (fTTj) comes to 0) the magnetic susceptibility becomes singular and the following increase 
of q results in a change of sign for the susceptibility from negative to positive leading to the 
stability of the paramagnetic state. 

Fig. EH shows also the Ni magnetic moment as a function of the wave vector q of spin 
spiral. As one can see, the magnitude of the goes down upon increase of q reaching m = 
at the critical value of the wave vector q c . 
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FIG. 7: Wave- vector dependent spin susceptibility of paramagnetic Ni having lattice parameter a = 
6.65 a.u. together with local Ni magnetic moment as a function of wave-vector characterizing non- 
collinear spiral magnetic structure (a). The wave- vector dependent unenhanced (b) and enhanced 
(c) magnetic susceptibilities for paramagnetic Ni calculated for different lattice parameters. 



As is shown in Fig. [7)d, a decrease of the lattice parameter leads to a decrease of the 
unenhanced susceptibility due to the broadening of the energy bands of the 3d-states. 
This results in a decrease of the critical wave vectors q c until a lattice parameter is reached 
for which q c = 0. For smaller lattice parameters the ground state of Ni is the PM state. 
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FIG. 8: Wave- vector dependent enhanced (a) and unenhanced (b) spin susceptibilities of param- 
agnetic Cr. For comparison, the unenhanced susceptibility is plotted also at the panel (a). 

3. Pure antiferromagnetic bee Cr 



Results for the non-magnetic state of Cr having the AFM structure as a ground state are 
shown in Fig. |SJ Note that the antiferromagnetic order of Cr on the one side is a result of 
nearly-half filling of the d-band 28| (similar to Mn), that should result in a commensurate 
AFM structure. However, Cr exhibits also an instability with respect to an incommensurate 
spin-density wave (SDW) with the wave vector q ^(0, 0, 0.95), which is a result of the Fermi 
surface nesting. This leads to a singularity of the magnetic susceptibility at q w ^(0, 0, 0.95) 
of paramagnetic Cr. This SDW instability in Cr and the corresponding behavior of the 
momentum dependent magnetic susceptibility was discussed in the literature by several 
authors |30l-l32]. 

Our present results demonstrate that the calculation of the momentum-resolved magnetic 
susceptibility properly reproduce its q dependent features for Cr. The calculations have been 
performed for a lattice parameter a = 5.4 a.u. which is slightly smaller than the experimental 



16 



one (a ?a 5.44 a.u.). At this lattice parameter the PM state was found to be more stable than 
the AFM state. This allows us to observe the behavior of x(o) due to the Fermi surface 
nesting avoiding the influence of other singularities connected to the instability around 
q = n/a(0, 0, 1) with respect to the AFM state. 

Fig. [Sb shows a monotonous increase of the unenhanced susceptibilities with increasing 
wave vector q reaching its maximum at q fa 7r/a(0, 0, 0.92). The enhanced susceptibility, 
also increasing with wave vector q, has a drastic increase at q ~ 7r/a(0, 0, 0.92) due the 
enhancement factor (Eq. (fl7|) ). which is associated with a singularity caused by the Fermi 
surface nesting mentioned above. 

Here, we do not discuss the q dependence of the exchange integral I(q) as this was done 
in detail by Sandratskii and Kiibler. Nevertheless, we would like to stress that this feature 
is taken into account within the self-consistent calculations for every wave vector. In fact 
this is essential for the analysis of the stability of the paramagnetic state. 

IV. CONCLUSION 

A theoretical approach for electronic structure calculations on systems with spiral mag- 
netic structures within the KKR Green's function formalism has been presented. As has 
been demonstrated, by making use of symmetry, the scattering path operator can be ob- 
tained by solving the corresponding equation of motion in the reciprocal space. Compared 
to the case of collinear magnetic structure only the structural Green's function to be used 
involves the wave vector of the spin spiral. As the KKR-formalism combined with the CPA 
allows to deal with chemically disordered materials, corresponding spin spiral investigations 
on various disordered alloys could be performed. In particular the energy of spin spirals and 
the behavior of the magnetic moments of the alloy components was analyzed. In addition 
is was shown that the approach presented can be efficiently used for the calculation of the 
momentum resolved longitudinal magnetic susceptibilities of pure materials as well as of 
disordered alloys. 
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